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By following the trajectories of quantum particles inside a periodic lattice and preserving their 
classical probabilities for reflection, transmission and absorption at each lattice plane, classical 
scattering outcomes are obtained. Diffraction phenomena in crystals are reproducible after assigning 
probability amplitudes to every classical outcome. When applied to X-ray diffraction in reflection 
geometry, this procedure has provided simple recursive equations to calculated the X-ray reflectivity 
of crystals with thickness varying since a few atomic layers to infinity. The results are in agreement 
with the dynamical theory of X-ray diffraction, even when absorption is considered. 
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I. INTRODUCTION 

Diffraction is essentially an undulatory phenomenon 
and, consequently, diffraction theories were developed us- 
ing a formalism that describes the propagation of waves 
through a periodic medium. The propagation of vec- 
torial electromagnetic waves is governed by Maxwell's 
equations while the propagation of wave functions is gov- 
erned by Schrodinger's equation. However, both type of 
waves are related to the amplitude of probabilities for 
possible trajectories of photons or matter particles (elec- 
trons, neutrons, ...etc) inside the medium. Therefore, by 
calculating the probabilities of all possible trajectories 
for classical particles, i.e. by visualizing the particles as 
globules, and assigning to each classical probability quan- 
tum probability amplitudes we must be able to achieve a 
complete quantum description of diffraction phenomena. 

This concept of equivalency between classical and 
quantum probabilities is further developed in this paper 
and applied to diffraction process in reflection geometry. 



II. GENERAL FORMALISM 

The outcome of elastic scattering of classical particles, 
like bullets, by partially permeable walls can be calcu- 
lated by introducing reflection and transmission proba- 
bilities. For example, two identical machine guns, sym- 
metrically displaced on each side of a wall, as shown in 
Fig. 1, are shooting at the frequency Nq. The numbers of 
bullets per unit of time reaching the Di and D2 counters 
are Ni = (r + t)No and N% = (f+t)No, respectively. r(f) 
and t(t) are reflection and transmission probabilities for 
bullets from the Si(S2) sources in which Ni + N2 = 2Nq 
since r + 1 = 1 and f + i = 1 . 

In the case of identical sources of quantum particles, 
ipR = Ripot ipT — Tipo, tpx — Rtpa and ipf = Tipo are the 
probability amplitudes assigned to each classical outcome 
given by the r, f, t, and t probabilities in Fig. 1. Since i/jq 
stands for the amplitude of the incident beams \iPr\ 2 + 




FIG. 1: Elastic scattering of particles from identical beam 
sources Si and S2, symmetrically displaced at both side of a 
partially-permeable wall. Classical particles only requires re- 
flection (r, f) and transmission (t, t) probabilities for predict- 
ing the outcomes. On the other hand, probability amplitude 
coefficients, R, R, T, and T are required for quantum particles 
such as photons and electrons; their square modulii provide 
the classical probabilities while their phase relationships obey 
Eq. 

hM 2 = fel 2 + h/M 2 = hAol 2 , or 

|i?| 2 + |T| 2 = |i?| 2 + |T| 2 = l. (1) 

As in the classical case, the sum of the beam powers 
(particle energy x frequency) at both detectors are iden- 
tical to the total power of the sources. It means that 
\ipR + V'tI 2 + + V"r| 2 = 2|-0o I 2 , and hence 

Rf* + R*f + RT* + R*T = 0. (2) 

This equation stipulates some phase relationships be- 
tween the reflected and transmitted amplitudes, as usu- 
ally obtained for laser beam splittersii Without loosing 
generality, the probability amplitude coefficients can be 
written as 



R = ±iVfe lC5+v \ and T = y/f (3) 
where according to Eq. 

|i?||f , |sin5+|i?||T|sin ( 5 = 0. (4) 

if and Cp stand for phase variation across the thickness of 
the wall] their values are arbitrary regarding convenient 
choices of references. The 8 and 8 phases are determined 
by the internal structure of the wall, which is in fact a 
plane of matter with thickness d, comparable to the par- 
ticle's wavelength A. When \R\ = \R\, and consequently 
\T\ = |f |, S = -6 according to Eq. gj. H 5 = 5 = 0, 
the reflected waves are 90° shifted with respect to the 
transmitted ones. 

Although the phase shift of reflected waves is a well- 
known phenomenon,!* 2 * we have reproduced its demon- 
stration here to emphasize the correlation between clas- 
sical and quantum probabilities in which the square mod- 
ulii of the probability amplitude coefficients are related 
to the classical probabilities as given in Eq. @ where 



ITT 



t, \R\ 2 = f, and ITI 2 = t. 



This concept can also be applied on a rather complex 
situation configured by one incident beam of particles 
and two parallel planes, as depicted in Fig. 2. Since the 
particles can suffer several bounces (reflections) at both 
planes before leaving the interplane region, classical re- 
flection and transmission probabilities are given by 



P, 



class. 



= r + trt + trrrt + trrrrrt + 



and 



pdass. _ t 2 + trft + trfrft + 



(5a) 



(5b) 



respectively. They correspond to the sum of probabilities 
for all possible outcomes of the scattering, in the sense 
that the particles can be scattered after one, two, or 
n-bounces. Each outcome can be identified as a classical 
n-bounce channel, where channels with odd (even) num- 
ber of bounces contribute only to the reflection (trans- 
mission) probability. 

The quantum coefficients of scattering are obtained by 
assigning to each classical channel, identified in Eqs. 
a probability amplitude with its respective phase. The 
phase depends on the total length of the trajectory of a 
given channel. If ip and if are chosen to be zero in Eq. 
the phase difference between two consecutive channels 
scattering on a same direction is an explicit function ei- 
ther of the plane thickness d, as well as of the separation 
between them h. However, if we choose 




FIG. 2: Classical n-bounce channels for the elastic scattering 
of particles in a two-plane system, r, f and t, i are classical 
probabilities of reflection and transmission for incidence on 
two different directions: from the top-left or from the bottom- 
left, e.g. sources Si and S2 in Fig. 1. [r, f] + [t, i] = [1, 1]. 



the amplitude coefficients will account for the phase dif- 
ference owing to the plane thickness. For instance, the 
probability amplitudes of channels 1' and 1 in Fig. 2 are 



i\)y = Ri/jq = ±iy/r e 



Wo 



and 



V>i = TRT1P0 = ±iVtrt e w e 



Jtp Md+<p) iip iAtp(h) 



^0- 



Then, 2ip + Acp(h) is the phase difference between two 
consecutive channels where Aip(h) = — 4irh$in6/A. The 
scattering time delay of higher order channels to with re- 
spect to the lower order ones is responsible for the minus 
sign of this expression, as well as in Eq. JHJ| . 

Of particular interest is the case where h — > 0. In 
such a case, the amplitude coefficients of two-plane sys- 
tems are obtained by direct replacements: r, f, t, t — > 
R, R, T, f in Eq. ©. It provides 



and 



Ro = R[l + TT(1 + RR + RRRR + ...)] 



T = T 2 (l + RR+ RRRR 



for incidence from the top-left, as shown in Fig. 2. In the 
other symmetrical situation, which corresponds to inci- 
dence from the bottom-left, e.g. source S2 in Fig. 1, the 
Ro and Tq coefficients for the reflected and transmitted 
amplitudes are obtained by analogous procedure. Since 



1/(1 -z) 



ii\z\ < 1, the probability amplitude coefficients can be 
written in a more compact format, 



^0 
Ro 



TT 



1-RR 



(7a) 



and 
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The i?o, T , R , and T coefficients also fulfil the 
energy /probability conservation equations, Eqs. and 
(J2J). Therefore, they provide a building block of recursive 
formulas for calculating the reflection and transmission 
amplitude coefficients for N-plane systems, like atomic 
planes in perfect crystals. 

In the present view, each individual scattering plane, 
with R, R, T, T coefficients, stands for an element of 
periodicity composed, in general, of several atomic planes 
that repeat themselves with period d in order to build 
up the crystal, as illustrated in Fig. 3. A fast recursive 
formula to go from very thin to semi-infinity crystals (of 
infinity thickness) is obtained by building the crystals in 
geometrical progression of N= 2™ planes (n = 1,2, ... ), 
and hence 
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FIG. 3: Elements of periodicity in crystals are composed, in 
general, of several atomic planes that repeat themselves with 
period d. The scattering properties of each element is given 
by reflection R, R and transmission T, T coefficients, Eq. . 
The intensity reflectivity of an array with N= 2™ repetitions 
can be calculated by Eqs. J5J. 



Rn 
R» 



and 



fn 



T n -iT n -i 

1 — Rn-lRn-l 



1 — R n -iR n -i 



Rn-1 
Rn-1 



T 2 

f. 2 



(8a) 



(8b) 



are the amplitude coefficients for a crystal of thickness 
t = Nd. It is obtained by analogy with the coefficients 
of the two-plane system in Eqs. Q. 

The maximum of the intensity reflectivity \R n (0)\ 2 , oc- 
curs around the diffraction condition of the lattice, where 
if = —mit and m is an integer number. The center of the 
first order reflection, m = 1, is observed at the Bragg 
angle 9b- 



A. Absorption 



t — a. Consequently, absorption is taken to reduce the 
transmission probability only, i.e. t = 1 — r — a in Eq. 

In this method, all diffraction channels of the two-plane 
system in Fig. 2, except channels 1', have the same ab- 
sorption probability of 2a. It corresponds to the intensity 
reduction ratio of a per interplane distance, and hence to 
an effective linear absorption coefficient fi — a sin 8 /d. It 
certainly is a valid approach when the interplane space 
is empty; for instance, when the elements of periodicity 
are made of single atomic planes. 



2. Method II 

When the interplane space is filled up with matter of 
linear absorption coefficient fj,, the absorption of a single 
plane can still be represented by the absorption proba- 
bility 



There are two methods to properly allow absorption 
in the medium. They differ on the distribution of the 
interplane matter - the matter that fills up the space 
between the ideal planes used to represent the periodicity 
of the structure. 



Method I 



\idj sin 9. 



(9) 



But, it does not account for absorption along of the zigzag 
routes of the channels (Fig. 2). Every time the particle 
crosses the interplane distance, its phase is delayed by ip 
and its exiting probability (from the interplane region) 
decreased by exp(— [idj sin 9). Then, the replacements of 
tp and (p in Eq. by 



Depending on the interaction between the particles and 
the atoms in the crystal, some of the particles will not be 
able to further participate on the diffraction process after 
reaching an atomic plane. It gives rise to general absorp- 
tion probability, a, accounting for the particles that are 
not reflected neither transmitted under diffraction con- 
dition, so that t + r = 1 — a. 

Since r and a are very small quantities compared to 
t ~ 1, the reduction in the reflection probability r(l — 
a) ~ r, due to absorption is negligible in comparison to 
the reduction in the transmission probability t(l — a) ~ 



/ 27r , • n ■ d 
(p = rising + i 

A 2 sin6» 



(10) 



provide the proper absorption in crystals. 

Since the ideal planes are responsible for the elastic 
scattering process only, the reflection and transmission 
probabilities are preserved in this method, i.e. t + r = 1 
in Eq. ©. Absorption is taken to reduce the square 
modulii of the coefficients so that 



\R[ 



re 



and IT! 



te 



t- 
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where e~ a ~ 1 — a. The same occurs to R and T. These 
relationships demonstrate that, the elastic scattering of 
each element of periodicity is not affected by absorption 
and, although, methods I and II seems to be physically 
different, they are in fact equivalent and provide the same 
results for a -C 1. Both methods can be used, depending 
on particular conveniences, but they must not be mixed. 
In more specific words, for a same absorption process, 
one must use t = 1 — r — a and <p [Eq. ©], or t = 1 — r 
and if' [Eq. QDJ]. 

A few examples on the behavior of the intensity re- 
flectivity, \R n (9)\ 2 , as a function of incidence angle, crys- 
tal thickness, reflection and absorption probabilities are 
shown in Fig. 4. No difference has been observed when 
using methods I or II for absorption. The maximum re- 
flectivity is always less or equal than 1 even for semi- 
infinity crystals. The transmitted intensities are given by 
|T n (0)| 2 , which is quite different than [I - |i?„(6>)| 2 ]e- QT 
in the center of the diffraction peak, as can be seen in 
Fig. 4 (top-left inset) . The phase of R n (8) is also shown; 
8 = 5 = in all cases. The analysis of the reflectivity 
profiles for thick low-absorbing crystals leads to an em- 
pirical equation for the full width of the half maximum 
(FWHM), 




W=-^/\R\\R\tan0 B . 



(11) 
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FIG. 4: Intensity reflectivity \R n (9)\ 2 , calculated according 
Eqs. © in crystals of N= 2" planes. A9 = 9-9 B , A = 1.54A, 
and d = 3.14A in all of the following cases: (1) [N, r, (i(cm -1 )] 
= [2048, 16 x 10" s , 0.01] (pink curve); (2) [4096, 16 x 10 -8 , 
0.01] (blue curve); (3) [oo, 16xl0~ 8 , 0.01] (red curve); (4) [oo, 
64 x 10" 8 , 0.01] ( green curve); (5) [oo, 64 x 10 8 , 200] (orange 
curve); (6) [oo, 64 x 10" 8 , 8000] (black curve); and (7) [4096, 
64xl0" 8 ,400] (dashed curve, top-left inset). \T n (8)\ 2 (closed- 
black circles) and [1 — \R n (9)\' 2 ]e'~ aT (open-gray circles) are 
compared in the top- left inset for case (7). For cases (3), (4), 
(5) and (6), the phases of R n (9) are shown in the inset at the 
top-right. 



B. Single bounce crystals 

As a final example of classical-quantum concepts of 
probabilities applied to diffraction physics, let consider 
a special crystal where back-reflections are not allowed, 
i.e. where f = and so t = 1 — a (absorption accord- 
ing method I). In this crystal diffraction occurs by inter- 
ference of probability amplitudes for a single scattering 
event in one of N planes. 

The classical probability is provided by 



N-l 



r> class. 
r sb 



r + trt + ttrtt + 



n=0 



1 - (tt) N 
' 1 - («) 



for t + r = 1 — a. However, in this case, for applying the 
same rule used before to convert a classical problem into 
its quantum counterpart, as from Eqs. J5J to Eqs. J7J); 
the probability amplitude of the incoming particle has to 
be extended to the entire lattice. 

To better visualize this, imagine the following situa- 
tion, N position detectors of area A n distributed, with- 
out overlapping, over a screen of area As- If a quantum 
particle has the same probability to hit any position of 
the screen, the probability amplitude at the ni/i-detector 
is ip n = \J A n jA s so that 



N-l 

E 

n=0 



1 N_1 



A-r, 



is the effective probability of measuring the particle in 
one of the detectors. 

By analogy, the amplitude reflection coefficient, R s b, 
for a single bounce in one of N planes has to be nor- 
malized by 



Ns, where N$ = y^Vj~/V| is the potential 
number of planes that the particle can reach on its way 
down and up inside the crystal. The reflection coefficient 
is then obtained from P s ^ lass - as 



R s 



r i - (Try 
TNs i-tt 



(12) 



R, T and T are given in Eq. (jSJ. For transmission prob- 
abilities t and t, the particle can reach a maximum of 
Ni = 1/(1 — t) and N-\ = 1/(1 — t) planes, respectively. 
With this normalization, the maximum possible value of 
|.R s b| 2 is 1, for < a <C 1, independently of the real 
number of N planes in the crystal as shown in Fig. 5. 



III. DISCUSSION ON X-RAY DIFFRACTION 

Diffraction of X-rays in crystals is the most common 
technique for studying how atoms and molecules organize 
themselves in the solid matter. To apply or compare the 
intensity reflectivity obtained from Eqs. (JSJ with X-ray 
results, scattering coefficients for the amplitude of the 
electric field are required. 
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FIG. 5: Intensity reflectivity \R s b\ , in single bounce crystals 
according to Eq. JTTJ- AO = 9 - 6 B , X = 1.5A, d = l.OA, 
r = 0.01, a = 0.002, and N = 60 (gray curve) or N = 10 e 
(black curve). The phase of R s b is shown in the inset. 



The scattering and absorption of X-rays by individual 
atoms can be described by 



It is possible to generalize the scattered amplitude of 
the single layer E$l, to elements of periodicity composed 
of several atomic planes by using the structure factor 
Fh, of the chosen reflection H. In symmetrical reflection 
geometry, the diffraction vector H, of reflection H, is 
aligned to the crystal surface normal direction, and hence 

M<r.(20) - ; 

V c 

V c is the unit cell volume and d = 1/\H\ as usual. With 
this substitution in Eq. (|14fl . the reflection coefficients 
are obtained as 



.r e X\C\d 

1 vTsmT 



(16) 



where |C| = 1 or cos 29 for a or ir polarizations, respec- 
tively. Note that in this development Fh is accounting 
for the elastic scattering only. It means that, 



F h = £(/ + /On exp{2mH ■ r„ 



(17) 



and 



a s (29)=r e X[f(29) + f']e 



cr = 2r e Xf" 



-B(28) 



(13a) 



(13b) 



when the X-ray polarization is linear and perpendic- 
ular to the scattering plane. The former provides the 
amplitude of the atomic elastic scattering with correc- 
tions for atomic resonance /', and Debye temperature 
factor B(29). f(29) is the scattering factor through an 
angle 29 from the incident beam direction, and r e = 
2.818 x 10 _5 A is the classical electron radius. The atomic 
absorption cross section cr a , depends only on /" correc- 
tion of the atomic scattering factor for absorption. More 
details on /' and /" are available on the International 
Tables for Crystallography. 

A single layer of atoms, all of the same kind, with M 
atoms per unit area scatters an incident plane wave of 
amplitude Eq according to 



E. 



SL 



M 



smf 



a s {29)E = R SL E 



(14) 



(Warren, 1969). The probability to measure the specular 
scattering of a single photon by this layer of atoms is, 
therefore, r = |i?si| 2 while the transmission probability 
is given by t = 1 — r — a (absorption method I, §11. A. 1.) 
where 



M 

a = - — -a n 



is the absorption probability per atomic layer. 



(15) 



n runs over all atoms in the unit cell, and temperature 
factors are now implicit in the / and /' values. 

Within the approximation of uniform absorption by 
the interplane matter, the linear absorption coefficient 



2tt 

T 



(18) 



where T = r e X 2 /irV c , is obtained from the average ab- 
sorption probability per element of periodicity in which 



Ma a 



2r e Xd ■ 
V c 



. f n 



in Eq. (|15fl : then Eq © provides the above fi expression. 

Since R and R are known, it is possible to compare 
the intensity reflectivity from Eqs. JSJ with that expected 
from the dynamical theory of X-ray diffraction2i2i4i5i£iL& 
in symmetrical reflection geometryAifliii The width of 
the total reflection domain, characterized by top-hat 
shape of the reflection curves, as in Fig. 4, also known 
as the Darwin width is, according to Authier & Mal- 
grangeji^ 



W D = ^^f^ = 3 -W. (19) 

7T V c COS 9 TT 

By replacing Eq. ((Toll into Eq. ltTT|) . the ratio between 
Wd and the width W, of the \R n {9)\ 2 curve, is observed 
to be very close to unit. The phase behavior of R n {9), 
inset Fig. 4, is also very similar to that calculated by the 
dynamical theory^ although the center of the reflection 
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curve at 9b are slightly different when A is corrected for 
the effect of refraction. 

The intensity reflectivity calculated here provides simi- 
lar results to the dynamical theory. The developed recur- 
sive formula, Eqs. ©, is quite equivalent to that devel- 
oped by DarwinpLiS although the actual deduction has 
been obtained by a simple procedure. It is almost as 
simple as that of the Laue's kinematical (or geometrical) 
theory ) 6 i 14 i 15 which was derived by adding the amplitudes 
of the waves scattered by each plane, and taking into ac- 
count only the optical path differences among them; but 
neglecting the interaction of the propagating waves and 
matter. It is possible to derive the kinematical reflec- 
tion coefficient from R s t, , Eq. (|12() , when t = i = 1 and 
Ng = N, which means that all planes contribute evenly, 
independently of how deep they are inside the crystal. 



A. On the physical meaning of a complex atomic 
scattering factor 

The seminal work of Darwir&iS followed by Prinsii 
gave rise to the actual method in which X-ray photoab- 
sorption is accounted in crystals. It is responsible for the 
/" correction and to the widely used expression of the 
atomic scattering factor / = / + /' + if", as a complex 
number. 

When calculating the index of refraction for an aggre- 
gate of atoms, the imaginary part of / does have a very 
clear meaning, which is to reduce the amplitude of the 
traveling waves in the medium due to absorption, as in 
absorption method II (§11. A. 2.). However, what would 
be its meaning when used, for instance, to calculate the 
scattered amplitude by the layer of atoms in Eq. (|14(l ? 

Rsl — —ir e \Mf/ sin# would lead to 

r = \Rs L \ 2 = ( r -^) 2 [(f+f) 2 + (rn 

\ smO ) 

which implies that the reflection probability increases as 
the absorption cross section increases since the latter is 
proportional to /". This increase in r owning to ab- 
sorption has no physical meaning; although in practice, 
this increase is very small (/") 2 *C (/ + f') 2 , the sum of 
the reflection and transmission probabilities is preserved 
t + r = 1, and the overall result is very close to that 
obtained by absorption method II. 

As a consequence of calculating reflection coefficients 
only with / and /', the structure factors in Eqs. and 
(|17|l do not break Friedel's Law, i.e. \Fjj \ = \Fh \ even in 
non-centrosymmetric crystals. Moreover, in the present 
description of the X-ray diffraction, based on preserva- 
tion of the reflection, transmission and absorption prob- 
abilities of X-ray photons, there is still another problem 
regarding the index of refraction. The atomic absorption 
cross section a ai Eq. i|13b(l . is so small that the arrange- 
ment of the atoms inside the unit cell does not affect the 
absorption probability of the X-ray photons when cross- 




FIG. 6: Griding of standing waves during rocking curves. The 
element of periodicity for the chosen reflection H is composed 
by the A and B atomic planes. The nodes of the standing 
waves are (a) passing over plane A for reflection H, and (b) 
over plane B for reflection H. The arrows indicates the move- 
ment of the nodes as the crystal is rocked from A8 — —Wd/2 
(solid lines) to A9 = + W D /2 (dashed lines). 



ing the elements of periodicity. Therefore, the only con- 
tribution of /" to be considered in the index of refraction 
is the average one that provides the linear absorption co- 
efficient in Eq. I|18|l . It implies that the absorption prob- 
abilities along the incident and reflected beam directions 
have the same value, and hence \T\ = \T\. 

All of the above facts are in agreement with the beam 
splitter equation Eq. where |i?||T| = \R\\T\ and the 
structure factor phases have opposite signals S = —5. 
Nevertheless, the intensity reflectivities of the H and H 
reflections can be slightly different. This difference has 
its root on the imaginary part of /, which is the /" cor- 
rection for atomic photoabsorption, as well as on the for- 
mation of standing wavesi&iLi^ inside the crystal. 

In reflection geometry, the pattern of standing waves 
formed by the interference of the incident and reflected 
waves was proposed by Battermani2i2S and it is a well- 
known phenomenon today In the present context, 
the nodes and antinodes of the standing waves deter- 
mine the positions where the photons have, respectively, 
the highest and the lowest probabilities to be found 
along their classical path inside the crystal. The absorp- 
tion probability is, therefore, weighted by the pattern of 
standing waves. 

As the crystal is rocked through the reflection domain, 
the phase of the reflection coefficient varies from 180° + 5 
to 5, as shown in Fig. 4 for 6 = 0. This phase variation 
induces the pattern of nodes and antinodes to grid by half 
a lattice plane distance inside the crystal, and then, the 
photons have higher probability to be absorbed at only 
one half of the element of periodicity. For the H and H 
reflections, the one halves scanned by the nodes of the 
standing waves are not the same and, consequently, the 
effective absorption for the Friedel's pair can be different. 
This hypothesis is illustrated in Fig. 6. 

To estimate the effective absorption coefficient ^ e // 
during a rocking curve, the intensities of the standing 
waves as a function of depth are approximated to 
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FIG. 7: Intensity reflectivity accounting for absorption mod- 
ulation due to standing waves inside a Ge crystal. A = 1.54 A, 
/ + /' = 30.9, /" = 0.92, and fj, = 352cm" 1 . From 
A9 = — Wd/2 to +Wd/"2, the intensity curves are calculated 
by using /j, e ff(9), Eq. 1211 . instead of fi in Eq. 1101 . The 
effective absorption coefficient through the Darwin width is 
shown in the inset. 







l + cos(^ + *) 


2 





~ /o ( z)cos 2^ + |^ (20) 

where 9 is the phase of the reflection coefficient R n , ac- 
counting for the 180° phase shift (Fig. 4, inset) across 
the Darwin width plus the phase S of Fh- It has also 
been assumed that the crystal is thick enough to assure 
a comparable strength between the incident and reflected 
waves at all depth. Since I(z) is a smooth function with 
very small variation over the lattice period, 

Mf {9) = ^rJ2f> S 2 (nH ■ r n + 9/2). (21) 

n 

This equation is obtained from Eq. (|18fl when using the 
intensity variation of the standing waves, cos 2 (irH ■ r n + 
9/2), as a weight function for the atomic absorption cross 
sections at each 9 value. H ■ r n stands for the depth of 
the nth-atom along the interplane distance. 

During a rocking curve, the average effective absorp- 
tion due to modulation by standing waves is estimated 
by integrating fj, e ff(9) in the reflection domain, 

<Me//>=-/ Me//(*)d¥. 

By solving the integral we have that 

< Me// >= f - §r J^fn M^H ■ r n + 5). (22) 

n 



Note that the sum of < /i e // > and < /2 e // >, 
for reflection H, is equal to the linear absorption co- 
efficient fi. For instance, in a germanium crystal and 
A = 1.54A, n = 352cm" 1 while < fi e ff >= 255cm" 1 
and < p, e ff >= 97cm" 1 for the 111 and Til reflections, 
respectively. Since the average absorptions are not the 
same, the intensities of the H and H reflections can be 
different. 

However, the variation of fj, e ff(9) across the reflection 
domain also changes the profile of the intensity reflectiv- 
ity curves as demonstrated for a few examples in Fig. 7. A 
detailed review on X-ray absorption and standing waves 
in reflection geometry can be found elsewhere^ 

IV. CONCLUDING REMARKS 

The theoretical foundations of diffraction physics were 
settled since the beginning of the 20th century. At that 
time the quantum of light was just a concept proposed 
by Einstein to explain photoelectric effect, and it was not 
experimentally proved before 1925, in the same year that 
Heisenberg developed the basis of the quantum mechan- 
ics. Therefore, the Maxwell's electromagnetic theory was 
the only possible treatment for the X-ray diffraction in 
the early 1900s. Ever since, a huge number of articles 
and books have further developed and reported the seed- 
ing theories of X-ray diffraction as essentially an undu- 
latory phenomenon. Diffraction of matter particles like 
electrons and neutrons has received a similar treatment 
although by the Schrodinger's wave equation. 

In nowadays it is well accepted that in the helm of 
quantum reality probability amplitudes rule the existence 
of the particles of light or of matter. The most fundamen- 
tal distinction between classical and quantum phenom- 
ena are the interference of probability amplitudes that do 
not exist in the former. By following the trajectories of 
quantum particles inside a periodic lattice and preserv- 
ing their classical probabilities for reflection, transmis- 
sion and absorption at each element of periodicity, series 
of classical outcomes are obtained. Diffraction in crystals 
are reproducible after assigning, according to some previ- 
ously established criteria, probability amplitudes to every 
classical outcome. When applied to X-ray diffraction in 
reflection geometry, this procedure has provided simple 
recursive equations to calculated the X-ray reflectivity of 
crystals with thickness varying since a few atomic lay- 
ers to infinity. The results are in perfect agreement with 
the dynamical theory of X-ray diffraction, even when ab- 
sorption is considered. It offers an alternative description 
of the diffraction phenomenon, with a different point-of- 
view of the usual wave diffraction theories. 
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